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I. INTRODUCTION 

Entanglement is the key resource of the vast fields of 
Quantum Information Processing, Quantum Computa- 
tion, and Quantum Technology, for an introduction see 
e.g. dill]. For example, applications of entangled states 
are those for quantum key distribution [3], quantum 
dense coding and quantum teleportation [5]. Thus 
both, the identification and the quantification of entan- 
glement, play a mayor role for future applications [5]. 

The phenomenon entanglement is closely related to the 
superposition principle of quantum mechanics. A pure 
separable state is represented by a product of states for 
both systems. A general pure state is a superposition of 
factorizable states. The minimal number of such global 
superpositions denotes the Schmidt rank (Tj. A sepa- 
rable mixed quantum state is a convex combination of 
pure factorizable quantum states [7|. The generalization 
of the Schmidt rank to mixed quantum states delivers 
the Schmidt number (SN). This generalization and the 
introduction of SN witnesses is given in [5HTU]. The ex- 
perimental construction of states with a certain SN has 
been realized in The SN of a mixed quantum state 
fulfills the axioms of an entanglement measure, cf. [TH- 
[2] . More precisely, it is a convex roof measure as defined 
in [15j[T6|. 

The identification of entanglement in terms of entan- 
glement witnesses has been introduced in [T7]. For a 
given witness an optimization can be performed |18| . Re- 
cently, we proposed optimized, necessary and sufficient 
conditions for the detection of entanglement [19 . Note 
that the latter optimization and the optimization of an 
entanglement witness are inherently different. Our opti- 
mization procedure delivers so-called separability eigen- 
value equations. They resemble the well-known eigen- 
value problem, but they include the factorization prop- 
erty of pure local quantum states. 

In the present contribution we study the identification 
of the SN of a given quantum state. We provide a method 



that deliver all optimized SN witnesses. The optimiza- 
tion yields equations which will be discussed in two forms. 
These generalizations of the separability eigenvalue prob- 
lem delivers similar equations for arbitrary mixed SN 
states. Our method will be compared with the well 
known spectral decomposition and the Schmidt decom- 
position of quantum states. Properties of the solutions 
of these equations will be studied. With some funda- 
mental examples, we generate some general classes of SN 
witnesses. We apply them to mixed quantum states. 

The paper is structured as follows. We motivate our 
method in Sec. HT] In Sec. IIIII we reformulate the detec- 
tion of the SN of a quantum state by witnesses in terms 
of arbitrary Hermitian operators and optimized, neces- 
sary and sufHcient conditions. The optimization proce- 
dure will be discussed in Sec. |IV[ where we express the 
optimization problem in terms of a perturbed eigenvalue 
problem and discuss properties of these equations. In 
Sec. |V] we solve these equations for a wide class of op- 
erators, including all one dimensional projectors and op- 
erators defined in continuous variables. We apply the 
method in Sec. |VI| to identify a SN greater than one and 
two, for the case of a phase diffused two-mode squeezed- 
vacuum state. A summary and some conclusions are 
given in Sec. |VII| 



II. MOTIVATION 

Let us consider the following experimental situation, 
cf. Fig. [T] We have a beam splitter with squeezed- 
vacuum states in both inputs. These states have the 
same amount of squeezing, but in orthogonal quadra- 
tures. The output is the two-mode squeezed-vacuum 
state \q), [q = ee**^ and < e < 1), 



|g) = v/r^^g'=|fc,fc). 



(1) 
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One output is disturbed in terms of a phase randomiza- 
tion, which equally randomizes the phase Lp between —5lp 
and +5ip. The scenario under consideration could be used 
for transferring one part of an entangled state through a 
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noisy channel. The sender keeps the other part of the 
state, e.g. in a delay line, such as an optical fiber. 
The measured state is given by 

2'5</? J-5v 



1 



+8^ 



5<p 



dip\ee"^){ee"^\ 



(2) 



with the local unitary operation C/((p) — X^feLo ^*'^'^l^)(^l■ 
In general, this state is a mixed quantum state in contin- 
uous variable systems. Due to the linear independence 
of the two-mode squeezed-vacuum states, the rank of the 
operator psg is, in general, also infinite. 
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FIG. 1: (color online) The experimental realization for the 
considered setup is the following. In both beam splitter in- 
puts, we have squeezed-vacuum states. The squeezing is con- 
sidered in orthogonal quadratures. The 50:50 beam splitter 
delivers an entangled two-mode squeezed-vacuum state. The 
detectors Da, D_b resemble the measurement of a, so far un- 
known, test operator L. The fluctuation of the distance be- 
tween the beam splitter output and the detector Da yield the 
phase randomization ±5^. The distance between the beam 
splitter output and the detector Ds is fixed. 



The two-mode squeezed-vacuum state has been con- 
sidered, for example, as a recource for quantum telepor- 
tation [2ni m] , quantum dense coding [55] , and quantum 
memories [23] . It has been shown that this state violates 
a continuous variable Bell inequality [24]. Pertubations 
of the pure two-mode squeezed-vacuum state have been 
studied, such as phase and amplitude damping [25j . or 
noise due to the transmission in optical fibers [2S]. Here, 
we focus on the SN of the phase randomized two-mode 
squeezed- vacuum state ps^p. 

In the case of zero phase randomization, we have an in- 
finite SN for po = |e)(e|. This state includes an entangled 
qubit (fc = 0, 1), an entangled qutrit [k = 0, 1, 2), . . . , an 
entangled qudit (fc = 0, 1, . . . , d— 1), see Eq Q. The two- 
mode squeezed- vacuum state is a global superposition of 



infinitely many product states fc), which corresponds 
an infinite SN. In this case, entanglement determines all 
the correlation between system A and system B. 

Now let us consider the fully randomized state dip = tt. 



Ptt = 



T 2 °° r+Ti 

k,l=0 



{l~e')J2e''^\k,k){k,k\. 



(3) 



k=0 



This state is obviously separable. Thus, it does not in- 
clude any entangled qudits, and the SN is one. The sub- 
systems A and B are only classically correlated. 

Some questions arise automatically. Somewhere in be- 
tween the two extreme cases. Sip — and Sip = tt, all 
the qudits disappear one after another, and the quan- 
tum correlations between A and B vanish. We aim to 
answer questions such as: For which value Sip the state 
becomes separable?, or: Which phase randomization Sip 
still deliver a quantum state with, for example, an en- 
tangled qutrit? With other words, we want to quantify 
the entanglement of the state ps^ under a certain phase 
randomization by the SN. Let us divide the problem into 
several sub-problems, which we solve in this manuscript: 

a. Conditions. Wc will derive general, necessary 
and sufficient conditions for the identification of a cer- 
tain Schmidt number r (SN-r). This method applies to 
all mixed quantum states. These conditions are based on 
measuring a general observable L, (L). For obtaining the 
SN of the state, we relate (L) to the maximal expectation 
value under all quantum states with a SN less or equal 
to r, fr{L). 

b. Optimization. We derive equations in two differ- 
ent forms, to obtain the desired function f^.. These equa- 
tions resemble the eigenvalue problem. The global super- 
position property of SN-r quantum states is encoded in 
these equations. 

c. Solution. We analyze these new kinds of equa- 
tions. We find a close relation to the Schmidt decom- 
position and the spectral decomposition of an operator. 
Afterward, we solve them for a large class of operators - 
including those for our considered example. 

d. Applying the method. We apply our method to 
the phase randomized two-mode squeezed-vacuum state. 
We show how the squeezing and the phase randomization 
influence the quantum correlations between subsystems 
A and B. We identify entanglement in general and the 
existence of a qutrit in the randomized state. 

Altogether, we derive a new mathematical method for 
the identification of the SN. This method allows us to 
test up to which order entangled qudits are contained in 
an arbitrary bipartite quantum state. These entangled 
qudits can be used for quantum information processing. 
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III. SCHMIDT NUMBER STATES AND 
WITNESSES 

Let us consider a bipartite quantum system which is 
given by compound Hilbert space V. = Ha ^ Hb- Previ- 
ously we have shown in, how entanglement in cantinuous 
variable systems can be identified in finite dimensional 
spaces |27j . Thus, it is not a restriction if we assume 
dim "H = dAds < oo. In addition let us denote the sets of 
linear operators acting on the Hilbert space as L\n{H, H), 
and the Hermitian operators as Herm(?^). 

Statistical mixtures of pure factorizable states define 
the separable quantum states Following this idea, 
quantum states with a SN less or equal to r are mixtures 
of pure states with a Schmidt rank less or equal to r. 
Pure states l^') with a SN less or equal to r, r{ij;) < r, 
are elements of the set Sr^^'^'^^ and can be written as [1] 

\ij) = U(g>vY,^k\k,k)^J2^k\ekJk), (4) 

k=l k=l 

with a local unitary operation U ^ V, orthonormal |fc), 
|efc), and \fk), and Schmidt coefficients Afc > 0. The 
mixed SN-r states a are elements of Sr given as con- 
vex combinations of those pure ones, tS^^ — convtSr^ ^. 
Let us note that the maximal possible SN is d = 
min{(i^, d^}, and the minimal one is r = 1 for separable 
quantum states. 

Analogously to entanglement witnesses, a SN witnesses 
is a Hermitian operator, with 

ya eSr -.TtctW >0, (5) 
3geSd: Tr gW < 0. (6) 

An operator fulfilling Eq. ^ is called SN-r witness, W G 
yVr- The condition given in Eq. (|6| will be considered at 
the end of our treatment. Here, let us only note that the 
set Wr includes all SN-r witnesses fulfilling Eq. and 
all positive semi-definite operators. Here, a SN witness 
W is considered to be optimal, if there exits a SN-r state 
a eSr with Tr aW = 0. 



A. Conditions for Schmidt number states 

In this subsection we aim to derive the conditions for 
the identification of SN-r states. It has been considered 
in [51-fTU]. that a quantum state has a SN greater than r, 
if and only if 

3W eWr-. Tr gW < 0. (7) 

The existence of such a witness is a consequence of the 
Hahn-Banach Theorem. Whereas, the general structure 
of SN-r witnesses - elements of Wr - is not clear. Thus, 
it is an advantage to find conditions in terms of arbitrary 
Hermitian operators with a well-known structure. 



For this purpose let us define the following function fr- 
This function maps a Hermitian operator to its maximal 
expectation value of all mixed and pure SN-r states. 

Definition 1 The maximal SN-r expectation value is 
given by 

/,(L) = sup{(V'|i|V'):|V')(^|e5(P"-^}. 

Proof. This value exists due to the compactness of 
^^purc)^ This value is the maximal expectation value of 
all SN-r states: 

"<": cSr-^ fr{L) < sup{TrCTL ; ct G 5 J 

= E|^,)(^,|e5(p"-) G Sr, withp^ > and 

fr(L) ^^p^fr{L) 

i, 

fr(L) > sup{Tr crL -.aeSr} 

It follows fr{L) — sup{Tr cri : a G Sr}- B 

In the following, the decomposition of all operators of 
the set yVr in terms of arbitrary Hermitian operators is 
considered. We show that a certain form delivers a SN-r 
witness, and all SN-r witness can be written in this form. 
This yields the identification of the SN. 

Let us consider an arbitrary Hermitian operator L G 
Herm('H), and a number A > fr{L). We define the oper- 
ator Wl as 

Wl = XI(^l- L, (8) 

with the bipartite identity operator I® I. This is a gener- 
alization of the construction of entanglement witnesses, 
cf . [ini [2H] • Obviously holds for all SN-r quantum states 
a 

TraWL = XTia- TrcrL > A - /^(L) > 0. (9) 

Thus, Wl is a SN-r witness. 

In the case A — fr{L), we call this SN-r witness opti- 
mized. For r = 1, our optimality condition is, in general, 
weaker then the definition in Ref. [IS]. To be precise, our 
optimized witness, A = fr{L), is finer than the witness 
for A > friL). 

The other way around we decompose all witnesses in 
this form. Let W G Wr be a SN-r witness. Using Eq. ([s]) 
multiplied with —1 and the resulting choice A = 0, the 
operator L — — W delivers the decomposition, 

= A > fri~W) ^ W-w = - {~W) = W. (10) 

It follows that Wl is a SN-r witness, and any SN-r wit- 
ness can be written as given in Eq. ([8]). This means the 
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set Wr and the set which is defined by operators of the 
form in Eq. ([s]) are identical. 

Moreover, the decomposition in Eq. ([s]) and ^ proves 
us the fact that optimized witnesses are necessary and 
sufficient for the detection of a SN. Summing up all these 
findings, we can formulate optimized, necessary and suf- 
ficient conditions in terms of arbitrary operators. 

Theorem 1 A quantum state g has a SN greater than r, 
if and only if there exists L € Herm('H) with 

(L) =TTgL> fr{L). 

Proof. The first part of the proof has been done above 
with the general structure of a witness in Eq. Q . We use 
this general form of an optimized SN-r witness Wl — 
fr(L)I^l — L. In addition, we use the fact that a witness 
must exist: 

> Tr (jWl = /r(i)Tr g - Tr gL ^ Tr gL > fr{L). U 

This condition means that the expectation value of g ex- 
ceeds the boundary given by all SN-r states. The formu- 
lation of the theorem could also be given by 

Tr gL < inf{(^|L|V') : \^){^\ £ ^^p^-^}, (11) 

if we use the minimal value instead of the maximal value 
in Definition [T] 

In conclusion, we have obtained SN-r conditions in 
terms of arbitrary operators. So far, only some special ex- 
amples of SN-r witnesses have been known, so that a gen- 
eral approach was unknown. This shortcoming we have 
resolved in terms of the condition in Theorem [TJ Such a 
method is of some general interest, since the identifica- 
tion of the SN of a quantum state delivers the structure 
and the quantification of entanglement in the notion of 
global superpositions. In the following, we will focus on 
the determination of f^, which is desired for a general 
characterization of entanglement through the SN. 



IV. THE r-SCHMIDT EIGENVALUE PROBLEM 

Now, we consider the function /,. as given in Defini- 
tion [I] The value of denotes the maximal expec- 
tation value of a given operator L under all SN-r states. 
First, we consider all extrema - not only the global max- 
imum. Obviously the largest extrema is the desired value 
of fr- The included optimization problem will lead us to 
equations which are of an algebraic nature. This means 
we transform a optimization problem to an algebraic sys- 
tem of equations. 



A. The optimization problem 

First of all, let us consider a weaker decomposition 
than the Schmidt decomposition given in Eq. (|4|. A pure 



quantum state {i/j) has a SN less or equal to r, if it can 
be written as 

r 

IV) =5^|xfc,2/fe), (12) 

k=l 

with {tp\ip) = 1. Here, the product states \xk,yk) are 
not necessarily orthonormal or linear independent. Some 
of them can also be the zero. Thus, lip) might have a 
Schmidt rank less than r. But in any case it follows that 

|V)(V|e5(P";'=\^ 

The optimization problem given in Definition [T] is the 
following. For a given L G Herm('H), we want to find the 
maxima or minima - in general extrema - of {iplLlip). For 

a quantum state |V')(V'I G Sr^™'^^ the expectation value 
is given as 



G(|V))=G(|xi),...,|x,),|yi), 
= (V'|i|V) ^ Gopt, 



|y.)) 



(13) 



where Gopt denotes an optimal value. We also make use 
of the normalization condition 



c(|v)) Mi^m - 1 = W i|^> -1 = 0. 



(14) 



An optimization problem under a certain condition can 
be solved with the method of Lagrangian Multipliers g. 
This means for all k holds 

dG dC 

91^=0, (15) 



d{xk\ 
dG 



d{xk\ 
dC 



d{yk\ d{yk\ 



= 0. 



(16) 



Using in the form of Eq. (12), we obtain from 
Eqs. (1151) and (fiel) for all k 



J2{yk\L\xj,yj) =g^{yk\yj)\xj), (17) 

j 3 

^{xk\L\xj,yj) =g^{xk\x.j)\yj). (18) 

j 3 

The Lagrangian Multiplier g can be obtained by summing 
up Eq. ( 17 ) - summation given by {xk | . . . ~ and using 
the condition C in Eq. (fT4|: 



9 ■ 1 ^ 9^{xk,yk\xj,yj) 



= ^{xk,yk\L\xj,yj) = Gopt. (19) 

k,3 

This is already the algebraic conversion of the optimiza- 
tion problem. 

Further on, let us simplify and rewrite these equations. 
Therefore we may consider the following definitions 



\y) 



(\yi) 



\ \Vr) 



Ls = (tlAL [\xk){Xj\ <^^)k.j = l,...,r ■ 

{tTBL[l<E>\yk){y,\])^^^^^_^. 



Lg 



(20) 

(21) 
(22) 



5 



Analogously to the operators Lg and L.j, we define the 
pseudo-metrics Ig and From Eqs. (17) and (181 we 
obtain the first set of equations for the SN^ optimization 
problem, 



Ly\x 



(23) 
(24) 



This form generalizes the separability eigenvalue equa- 
tions, as derived in [19]. In fact, for r = 1 



Eqs. (23) and (24 1 reduce to the separability eigenvalue 



equations. 

Definition 2 The r-Schmidt eigenvalue problem is de- 
fined in Form 1 as 

Ly\x) ^ gly\x) and Ls\y) = gls\y). 

We use the abbreviation r-SE for r-Schmidt eigenvalue. 
The vector \ip) is the r-Schmidt eigenvector (r-SE vec- 
tor), and the value g is the r-SE for \^|;) of L. 

The conversion of the optimization problem to its alge- 
braic form is given in the following theorem which sum- 
marizes the above calculations. 

Theorem 2 The SN-r state \ip) delivers an extremal ex- 
pectation value g — {ijj\L\^j) , if and only if it solves the 
r-SE problem defined in Definition^ ■ 

More sophisticated is the derivation of the second form 
of these equations. We note that the SN-r state in its 
weak form of Eq. ( 12 ) must also have a strict Schmidt de- 
composition, see Eq. ([4|. Obviously the optimization and 
the optimal value Gopt — g — {'4^\L\ip) does not depend 
on the form of a decomposition of the pure state. Let us 



note that the r in Eq. ( 12 ) is not necessarily the Schmidt 
rank r > r{'il)) of | 



However, due to the independence 

.(pure 



of the decomposition and !■(/;) (V^j G 



Hence, we relabel formally 

\xk) |efc), \xk) Afcl/fe) and r(V') ^ r. (25) 
Now, let us reconsider which information is included 



in Eqs. (17) and (18 1, and rewrite them as 



{ek\L\'4>) =5(efc|V') = 9>^k\fk), 
XkifklLW ^gX,{f,\^j) ^ gXl\ek) 



(26) 
(27) 



We use the notion {ekl^}) £ V-b and G for 

the projection of the state lip) on the \ek) and |/fc) 
component, respectively. Note that > 0. Thus, 



for each fc, we can multiply Eq. (27 1 with X^. . Equa- 



tions ([261 and ([27]) yield for k, k' . . . 

{ck, fk'\L\ip) = gXkSk,k'- 



(28) 



On the other hand, let us consider the action of L on 
the pure SN-r state {ip). This delivers the (not normal- 
ized) result IV''): 



with the decomposition of {ip') into a parallel component, 
5 IV'), and an orthogonal one, \x)- Now on the subsystems 
A and B act \ek) and \fk'), respectively, 



{ek,fk'\L\ip) ^{ekJk'W) 



We combine Eqs. (28) - (30) and obtain 



{&k,fk'W) =5(efc,/fc'|V^) + (efc,/fc'|x) ^9{ek,fk' 



(30) 



(31) 



for all k, k' = 1, . . . ,r. We conclude that \x) is orthogonal 
to each \ek,fk') for /c, fc' = 1, . . . , r which leads to 



fc— r+l k'—r-]-! 



\ek, fk' 



(32) 



This means that \x) is not only orthogonal to |V'), but 
also to all linear combinations of \ek, fk') £ Ha ^Hb oi 
the Schmidt decomposition of {ip) and k,k' — l,...,r. 



We denote this property of Eq. (32 1 as bi- orthogonality. 

From these considerations follows another algebraic 
conversion of the optimization problem as 



(33) 



LW^g\rP) + \x), 

with a bi-orthogonal vector \x) as it is given in Eq. ( [32| . 
We define the second algebraic form of the optimization 
problem for a state with maximally r global super- 
positions of factorizable states. 

Definition 3 The r-Schmidt eigenvalue problem is de- 
fined in Form 2 as: 

L|V)=5K') + lx), 

with a bi-orthogonal perturbation \x). 

Due to the fact that Definitions [2] and [3] solve the same 
optimization problem. Form 1 and Form 2 of the r-SE 
problem are equivalent. 

Theorem 3 The SN-r state \ip) delivers an extremal ex- 
pectation value g = {ip\L\^) , if and only if it solves the 
r-SE problem defined in Definition^ ■ 

B. The r-SE problem — Preliminary conclusions 

So far, we have derived two forms of the r-SE problem. 
Our initial intention was to determine the value of the 
function fr{L) for a given linear operator L. The r-SE g 
is the extremal value of the function G, 



G(|V^)) = (V'|L|V)=Gopt = g. 



(34) 



Under all extremal values the global maximum can be 
found as 



^IV^) = lV'')=5lV') + lx), 



(29) 



Corollary 1 fr{L) — max{g : g r-SE value of L}. 
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In some cases L cannot be used for the identifica- 
tion of SN-r states. Obviously holds that L is not suit- 
able in such a case, if and only if fr{L) = fd{L), cf. 
Eqs. Q and ([6| and d denoting the maximal possible 
Schmidt rank. This means, no quantum state with an 
arbitrary SN exceeds the maximal expectation value of 
L for SN-r states. Thus, the test operator L for the 
identification of the SN must have the following prop- 
erty. The eigenspace of the largest eigenvalue does not 
include a vector with a Schmidt rank less or equal to 
r. For r = 1 this is related to the range criterion for 
quantum states |29| . 

It is clear that the bi-orthogonal perturbation \x) has a 
Schmidt decomposition as well. Due to the bi-orthogonal 
property it can be given as 



(35) 



k=r+l 



Now the r-SE equation in Definition [3] can be rewritten 
as 

r r d 

L^Xr\ekJk)^g^K\ek,fk)+ ^ >^k\ek,fk)- 



fe=i 



fe=i 



k=r+l 



(36) 



Thus, we can write a necessary condition for a solution 
of the r-SE problem. 

Corollary 2 A r-SE vector \tp) and the vector Llip) have 
a Schmidt decomposition with the same local unitary op- 
eration U (>^V , 



k=l 
d 

id {U®V)L\i;) ^^\k\k,k). 



k=l 



This condition is only necessary. It would also be sufh- 
cient, if Aj, = gXk for fc = 1, . . . , r. But Corollary [i] in its 
given form is already a quite strong restriction for finding 
r-SE solutions. We will see this later when we apply our 
method. 



C. Relation to the eigenvalue problem 

We consider the r-SE problem for all possible r = 
1, . . . ,d. In addition, let us consider the (ordinary) eigen- 
value problem of L g Lin('H,H) in linear algebra. 



(37) 



with the eigenvector 1-0) and the eigenvalue g. It is clear 
that in this case there is a trivial perturbation \x) — 0. 
Thus, eigenvalues and eigenvectors, g and 1-0), are also 
solutions of the r-SE problem for r = r{'4>). The other 



way around, we conclude for r = d that no bi-orthogonal 
perturbation can remain. Therefore, we obtain in this 
case that the d-SE problem and the eigenvalue problem 
are identical. 

Now, let us consider local invertible operations. These 
are operations which can be written as S<^T, with S and 
T invertible operations, S'^S = I and T'^T = I. These 
operations cannot change the amount of entanglement, 
given by the SN, for any quantum states [301111]. Now 
we transform the r-SE problem of the operator L given 
in Definition [3] to the r-SE of an operator L' 



L' = {s (g) T) L (s (g> ry 

L'm=g'm-^\X'), 



(38) 
(39) 



This can be done by multiplying the r-SE problem of L 
with S* (g) T, 

{S ® T) = (5 (g) T) L (5* ® T)-^ {S (g) T) |0) 

^giS®T)\^) + {S®T)\x), (40) 

We obtain that the r-SE values remain, and the vectors 
are transformed by S' g) T, 

g = g', m = {S^T)\i,), and \x'} ^ {S ® T) \x} . 

(41) 

Corollary 3 Let L' be a locally transformed operator 
L' = (S <Si T) L (S T) , g is a r-SE value and |^) a is 
r-SE vector of L. It follows that 

g = g' and \^') ^{S®T) \^) 

are r-SE value and r-SE vector of L' , respectively. 

Proof. Any local invertible map can be decomposed 
in terms of unitary and diagonal matrices, e.g.. by the 
singular value decomposition. First let us consider the 
local unitary U (C/ g) V)~'^ = W (g) We use 

U ®V\ek, fi) ~\uk,vi) (orthonormal basis) 
L' ==([/g)y)L(C/1' g)0), 

with the r-SE problem 

r r d 

L^Xk\ek,fk)=g^Xk\ek,fk)+ ^ Afc|efe,/fc). 

k = l k=l k=r+l 



It follows from Eq. ( 40 1 that 



L'^Xk\uk,Vk) =g'^Xk\uk,Vk) + ^ Xk\uk,Vk). 

fc=l fc=l k=r+l 

In addition we consider local, diagonal, and invertible 
maps, 

Da<E) DslckJi) =^dk,Adi,B\ek, fi) {dk,A,dLB 7^0) 
L' ={Da g) Db)L{D^^ g) Dg'^), 
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and we conclude for the transformed r-SE problem 

r r 
L' ^ dk,Adk,B^k\Rki Ik) —g ^ dk^Adk,B^k\ekT fk) 



fc=l 



k=l 



dk.Adk,B^k\ek, fk)- 



k=r+l 

The Schmidt coefficients change, i-> jV''): but the r- 
SE value remains invariant. ■ 

For the ( ordinary) eigenvalue problem in linear algebra 
an invertible - in general global - transformation Tab can 
be found for the diagonalization of the matrix. In general, 
a matrix L can be transformed into another one by L' — 
Tj^sLT^g. This transformation changes the eigenvectors 
in the same way, and the eigenvalues remain invariant. 
For the r-SE problem we have a related situation which, 
in addition, reflects the global superposition property of 
SN r states. 

Let us consider a special case of local invertible opera- 
tions, cf. proof of Corollary [3j The same considerations 
are obviously true ioi S ®T — U with U and V be- 
ing unitary operations. Consequently, the r-SE problem 
is invariant under local unitary, which is of importance 
for the quantification of entanglement, cf. e.g. [2l|6]. 



V. SOLUTIONS OF THE r-SE PROBLEM 

So far we have considered general properties of the r- 
SE problem. We have compared our method with the 
eigenvalue problem in linear algebra. Now let us apply 
our method to some examples. Starting from low rank 
matrices, we explain how the r-SE problem can also be 
solved for higher rank operators. We consider an ap- 
proach connecting the spectral decomposition with the 
r-SE problem. Our SN-r conditions will be tested for 
some examples in continuous variable systems. 



A. One-dimensional projections 

For simplicity let us here and in the following restrict 
to Hilbert spaces with dim Ha = dim Hb = d. We 
consider a one dimensional projection L = with 
1 < r((^) < d, 



r(0) 



^Kk\k,k), 



Kk G C\{0}, and the normalization J2k I'^feP = 1- 
dition, let us assume that the coefficients Kk are ordered 
as |ki| > |k2| > • • • > l'^r(0)l- It is also useful to define 

Kr(^) + l = • • • = = 0. 

We aim to solve the r-SE problem for L. First, there 
are a lot of trivial solutions given by orthogonal eigen- 
vectors 

=0|V) + |x), for IV') ^10) and |x) = 0. (43) 



An example is lip) = (|2, 1) -I- |l,2))/\/2 for r = 2. 

The non-trivial solutions, g 7^ 0, are more interesting 
for obtaining fr{L). Applying Corollary [2] we conclude 
that we only need to consider states with a decom- 
position as 



IV) =^Afe|A:,A), 



(44) 



with Afe G C. The Schmidt rank of |V') is the number of 
coefhcients with 7^ 0. For a Schmidt rank less than d 
some of the \k have to be zero. 

Thus, we obtain by applying the projection, L\^j) = 
9\4') + Ix), 



L|V)=(</'|V)I0)= 5]<Aj^«:fe|fc,A) 

\l=l I k=\ 

d d 

^gY,Xk\Kk)+Y,\k\k,k). (45) 



k=l 



A closer look on Eq. ( 45 1 delivers all the possible solu- 
tions: 



A, 



CKfc or \k = 0, 



(46) 



with a constant number c — const. ^ 0. The bi- 
orthogonality property between and |x) delivers that 
from A/c 7^ follows A^ = 0. 

For the case r = 2 we obtain all solutions as 



IV'fe,;) 



1 



{Kk\k,k) +Ki\lJ))^ 



(47) 



gk.i = (ipkA^i'k.i) 



2\2 



{\^k\' + \^in 

\Kk\^ + \^l\^ 



\nk\ 



for k,l — 1, . . . , d and k ^ I. The maximal 2-SE vector 
and value are 



1 



1,2/ 



(kiIi,i) 



K2|2,2)), 



(48) 



51,2 = (MLW 



i\ni\' + \n2\Y 

klP + |K2p 



kip + k2p. 



In conclusion, every lip) - with r Schmidt coefficients 
Kfc = Afc 7^ - and the same Schmidt decomposition as 
(42) \<P) solves the r-SE problem. The r-SE value is given 
as g — |(0|'0)P- Using the Schmidt coefficients with the 
largest absolute value, we obtain fr{L) — J2k=i kfcP- 

For example, let us use r(0) = d — 00 and the projec- 
tion operator L constructed from 



= k) = yr^^g'^-|fc,fc). 



(49) 



fc=0 



for < l^l = e < 1 (note that the index starts from 0). 
This is the two-mode squeezed-vacuum state. In this case 
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we obtain the maximal solution of the r-SE problem as 

r-l 



9r = \{A\(t>)? 



1 



\2r 



(50) 

frmi^D- (51) 



For the quantmn state po = |e)(e| {q = e) we can identify 
the entanglement with an arbitrary SN by 



Trpo^ = 1 > 1 



(52) 



This optimal violation, Tr pqL = foo{L), does not depend 
on the amount of squeezing \q\. Even in the case of a 
minimal squeezing, \q\ <C 1, the condition Eq. (52) proves 
that the state po has a SN larger than any r. 

Another example could be r(0) ~ d and equally dis- 
tributed Schmidt coefficients, 



(53) 



fe=i 



In this case we obtain the maximal solution of the r-SE 
problem for L — as 



K) = ^^|fc,fc), 



(54) 



fe=i 



5, = |(V;,|0)|2 = -=/,(|0)(0|). (55) 

In conclusion, we can formulate the following necessary 
criteria. A quantum state g has a Schmidt rank greater 
than r if 



k=l 



(56) 



with \Kk\ {k — 1, . . . ,r) being the r largest Schmidt co- 
efficients of 16). 



Note that the singular value decomposition of resem- 
bles the Schmidt decomposition oi {ip) [IJ. This rewriting 
delivers that the scalar product of two states and |^) 
is the scalar product of the corresponding matrices 



(eiV) = Tr A^M^. 



(59) 



A general Hermitian operator L E Herm('H) can be de- 
composed as 

L =y^£fc|-0fc)(V'fc| 

= ^Lfc(Mfe®I)|$)($| (a4®i) , (60) 
fc 

with Lfc S M. An example would be the spectral decom- 
position of L with the eigenvalues Lk and the eigenvectors 
{Mk (g) I) |4>). We apply L on a given state (M^ I)]*) 
which yields 

L{M^^l)\<P) =Y,Lk (Tr [m^M^]) (Mfe(g)I)|$) 



LfeTr 



MlM^ 



Mk <»I|$). (61) 



From Corollary [2 follows that and the resulting 



matrix 



(Efc 



mIm, 



Mkj have the same singular 

value (Schmidt) decomposition, if (M^ (g)I)|$) is a r-SE 
vector. Note that the Schmidt rank of the state is 
the rank of the corresponding matrix, r{il)) = rankM^. 

Let us consider an example. May the matrices have 
a decomposition as M^. = UD}J/\ with unitary opera- 
tors U and V and a diagonal matrix with complex 
components. We have seen, that unitary operations ba- 
sically do not effect the r-SE problem, cf. Corollary |3] 
Hence, we choose U = 1 — V . The effect of this example 
is that the mapping delivers 



B. Higher rank operators 

An arbitrary Hermitian operator L can be decomposed 
in terms of one-dimensional projectors. In the following, 
we generalize our results for one-dimensional projections 
to classes of higher rank operators. It is of advantage to 
consider all one dimensional projections as a local trans- 
formation of a given state with maximal Schmidt rank. 
We start from pure vectors jV') and |<I>) 

d d d 

l^) = EE^'=.'l^'0> and = (57) 



fc=i 1=1 



The linear map — i^k,i\k){l\ yields the desired 
property 

IV') = (M^0I)|$) = (I® AfJ) 1$). (58) 



(62) 



which is already given as a diagonal matrix. In such a 
case it follows, that the vector L{M^ ® 1)1^) is given in 
Schmidt decomposition and we can treat this case analo- 
gously to those given in Section [V A| for low rank opera- 
tors. We may choose a decomposition Dk = X]j ^k,j\j){i\ 
(with Kfcj- £ C and |k| ^ \ — 1), and the operator L as 

fc 

= ^ LkKk,mKl^Jm,m){n,n\ 

k,m,n 

= ^lrn,n\m,m){n,n\. (63) 
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The new parameters 7m. n are given by 



fc 



(64) 



It is worth to note that this operator L is symmetric with 
respect to exchanging the quantum systems. 

For obtaining the nontrivial solutions, g ^ 0, we con- 
sider It/i) = A„|n,n). This yields 



= X! ^^lrn,nK \m,m). 



(65) 



Due to the fact that some can be zero, we obtain the 
r-SE solutions by neglecting some rows and the corre- 
sponding columns, cf. Eq. (46 1. These rows are given by 



the Schmidt coefficients oi[4>) which are zero, \q = 0. 
This means we restrict to the rows and columns with 
the index gi, . . . ,9^, and solve the (ordinary) eigenvalue 
problem of the resulting operator. 



i(gi,...,(j,r) = ^ iqi,q,\qi,qi){qj,qj\- (66) 

The simplest case is ?- = 1. We have to stroke out 
d— 1 rows and the corresponding columns. It immediately 
yields 



^(91) = iqi,qi\qi,qi){qi-,q2\- 



(67) 



The maximal expectation value for all possible choices of 
qi = n = 1, . . . , d delivers the function fi{L) as 



/i(L) = max7„^„. 



(68) 



In the case r = 2, we have to consider all principal 2x2 
sub-matrices of 7m^„ with (ji, 92 — ™, n. We calculate the 
maximal eigenvalue, and obtain the function f2{L) as 



/2(L) = max 



(69) 



m,7n ^n.n)'^ ~l~4:|7m,n| 



squeezed-vacuum state in Sec. V A This operator was 



suitable for the detection that the state po, 5ip = 0, has 
an infinite SN. However, we already observed that a total 
randomization, S^p ~ tt, delivers a separable state Pt,, cf. 
Eq (§. 

Now, we apply our considerations to the general phase 
randomized (mixed) state pstp to answer our initial ques- 
tions. Namely, for which values of 5lp the entanglement 
survives (SN r > 1), and which values of 5ip guarantee 
that the state contains an entangled qutrit (SN r > 2)? 
A phase randomization in a single mode state has been 
experimentally realized, and the non-classicality of such 
a state has been verified [33J. Here, we consider a sim- 
ilar situation, for two modes, where the desired aspect 
of non-classicality is the number of global superpositions 
between these modes. The method presented in the fol- 
lowing also applies to a non-equally distributed phase 
randomization in both subsystems A and B, and an ad- 
ditional amplitude randomization of e, as it occurs in the 
turbulent atmosphere [33]. 

Here, we chose a test operator L of formally the same 
form as the state ps^p in Eq. ([2|, 



2 /•i5¥' 



dip q"^\in, m) q*" {n, n\ 

-^V m=0 n=0 



= lm,n\m,m){n,n 



m,n— 



7m,n =(1 - / d^e'^^^'-^K 

-Sip 



(70) 
(71) 



with q = ee"^ and < e < 1. For simplicity, we have 
assumed an equally distributed randomization for the 
angle, — tt < —dip < ip < 5ip < n. Therefore, the cor- 
responding parameters 7m,„ read as 



m+n 



2, sin {5ip[m — n]) 



(1 



(1-^ ) A [ 1 ■ 

Oip[m — n\ 

e2)e™+"sinc {Sp[m ~ n]) 



(72) 



First, we may consider a fully randomized and equally 
distributed phase for the test operator, 6ip = tt. 



denoting the maximal eigenvalue of all principal 2x2 
sub-matrices. Analogously, the method can be applied 
to SN-r states by finding the maximal eigenvalue of all 
r X r principal sub-matrices of ^rn.,n- 



VI. PHASE RANDOMIZED TWO-MODE 
SQUEEZED- VACUUM 

Let us come back to our initially considered state psip, 
cf. Eq. ([2|. This state was generated by a local phase 
randomization of a two-mode squeezed- vacuum state, cf. 
Fig. [1] For witnessing the SN, we already considered 
the example of a projection defined by a pure two-mode 



L 



Ed 



e2)e2"| 



TO, m){m, m\ 



(73) 



The spectral decomposition of L is given in factorizable 
states. Thus, fi{L) — fao{L), and the phase insensitive 
operator cannot detect any entanglement. We already 
observed this result for the corresponding state p.„ which 
is separable. 

Now we consider a given partially phase randomized 
operator, L. We obtain from Eq. (68) for separable states 
(SN r = 1) a maximal expectation value of 



/i(L) = max 7„,,„ 

m— 0,. . . ,00 



- maxe 



(74) 
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The maximal expectation value for two-qubit states (SN 
r = 2), of. Eq. yields 



ML)= max i-^f e^^ + e'- 



(75) 



using n = m + k and fc > 1 we obtain 
1 - e 



/2(i) 



max e^™ max f 1 + e^'' 



2 meN fc>l 



[1 _ £2fe]2 _|_ 4£2fegjjj(,2 ^^j^j^ 



(76) 



The maximum over m yields m = 0. The maximum over 
k has to be calculated numerically. Note that fi{L) < 
f2{L), and therefore L can detect entanglement. 

Now let us consider the state ps^p. The expectation 
value {L)s^ = Trp^^L of the observable L for the state 
ps,^ reads as 

oo 

(X),^ =(l-e2)2 e^^^+^^sinc^iS^im-n]) 

m,n— 

oo / oo \ 

=(1 - 6^)2 ^ e'" 1 + 2 ^ e^'^sinc^S^k) 



k=l 



l + 2^e2'=sinc2((5<y9fc) 



(77) 



k=l 



Again, this can be calculated numerically. We consider 
the case e = 1/3, which corresponds to a moderate 
squeezing of about 3 dB [32] . The entanglement condition 
for the choice of an operator L = ps^p can be formulated 

as 



Trps^L-fr{L) = {L)s^ - fr{L) > 0. 



(78) 



This means a positive number is a verification of a SN 
greater than r. 



1.00 



0.75 



0.50 



0.00 



-0.25 




180° 



FIG. 2: (color online) Here, the detected entanglement for 
a two-mode squeezed-vacuum state is given. We choose a 
squeezing of 3dB. In addition the state is manipulated by an 
equally distributed phase randomization with the interval of 
the length 2S(fi. The functions are normalized to the maximal 
value. The entanglement can be verified until Sip — 79° (blue 
solid function). We also identified a greater Schmidt rank 
than 2 for phase diffusion of Sif = 25° (red dashed function). 



178° and a SN greater than 2 close to 102°, cf. Fig. [3] 
Here, the chosen test operator L is given as 



L = 



100 

E 

m.n— 



sine {Sip[rn ~ n])\m,rn){n,n\. (79) 



This means that the state pg^ is entangled even for a 
phase randomization up to to the total phase randomiza- 
tion. The other interesting point is that the state surely 
contains a qutrit (SN greater than 2) for a randomiza- 
tion of more then ±90° in the channel of the receiver. We 
note that for such a phase randomization, the squeezing 
of the single-mode squeezed-vacuum state would vanish. 



VII. SUMMARY AND CONCLUSIONS 



In Fig [2] the identification of entanglement for angles 
of an equal distribution in phase for [—Stp, Sip] is given. 
The values are obtained by numerical calculation for the 
matrix components 7m, « with m,n < 100. This means 
that the given positivity, i.e. the entanglement, can only 
increase for indices > 100. Let us stress that we choose 
an arbitrary operator for the detection of entanglement 
and a SN greater than 2. There may exists operators 
which deliver entanglement also for randomization above 
the limitations of the chosen operator L. 

In addition, it is worth to note that a higher squeezing 
of the state delivers a higher sensitivity. For a realistic 
10 dB (e = 0.82) squeezing of the input states, cf. [55] . 
we can identify entanglement up to a phase diffusion of 



We have derived optimized, necessary and sufficient 
conditions for the detection of the Schmidt number of 
an arbitrary bipartite quantum state. These conditions 
have been formulated in terms of arbitrary Hermitian 
test operators as measurable conditions. We have shown 
that the optimization problem leads to the r-Schmidt- 
eigenvalue equations. We discussed the properties of the 
solutions and consequences of these equations in con- 
nection with entanglement and its quantification via the 
Schmidt number. For example, we have shown which op- 
erators can be used for the identification of the Schmidt 
number, and we have considered the relation to the eigen- 
value problem. We have solved these equations for a wide 
class of operators, namely for all one dimensional projec- 
tions and classes of higher rank operators in finite and 
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1.00 



0.75 



— TrpL-fl(L) 
■--TrpL-f2(L) 



S 0.25 



-0.25 
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120° 
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FIG. 3: (color online) Here, the detected entanglement is 
given for higher squeezing, 10 dB, and the test operator in 
Eq. (791. In the numerical calculation, the entanglement can 
be verified for a phase randomization up to 5ip = 178° (blue 
solid function). We also identified a greater SN than 2 for 
phase diffusion up to Stp = 102° (red dashed function). 



infinite dimensional systems. We have shown a direct 



identification of the Schmidt number for all pure states 
with our condition. 

We also applied our method to the identification of 
the Schmidt number of mixed quantum states, includ- 
ing those with an infinite rank and infinite Schmidt rank 
in continuous variable systems. Examples of test oper- 
ators are considered, which identify the Schmidt num- 
ber for a broad class of mixed states. To illustrate the 
method, we have studied the example of a phase diffused 
two-mode squeezed- vacuum state for a squeezing of 3 dB 
and 10 dB. We have identified entanglement (Schmidt 
number greater than one), and we could identify entan- 
gled pairs of qutrits in the state (Schmidt number greater 
than two) for different values of the phase randomization. 
Moreover, the influence of the amount of squeezing and 
the randomization has been considered. 
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